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Interplay of Superconductivity, Antiferromagnetism, and Pauli Depairing in CeCoIn 5 

Carley Paulsen*, Dai Aoki 1 , Georg Knebel 1 , and Jacques Flouquet 11, 

Institut Neel, MCBT Department, CNRS and Universite Joseph Fourier, BP 166, 38042 Grenoble Cedex 9, France 
1 SPSMS, UMR-E CEA / UJF-Grenoble 1, IN AC, Grenoble, F-38054, France 

Low temperature magnetization M measurements down to 70 mK have been performed 
through the superconducting upper critical field H C 2 of CeCoIns with the field H parallel to 
the c axis of the tetragonal crystal. The crossing through the upper critical field H C 2(Q) is 
marked by a sign change in the temperature dependence of M(T). A Fermi liquid law of the 
magnetization is observed in the superconducting phase while in the paramagnetic state above 
H C 2 a quasi-linear temperature dependence is obeyed. Fhe specific heat below H C 2 is therefor 
dominated by the increase of j(H) , the linear term of the specific heat, and is marked by 
a singularity at H C 2{0)- No indication of the so called Q-phase was observed below H C 2(in 
contrast to the case with H applied in the basal plane). 
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Investigations of the heavy fermion Ce-115 family has 
lead to spectacular discoveries of new phases with a 
strong interplay between magnetism and superconduc- 
tivity (SC) as shown in CeRhlns with the re-entrance 
of antiferromagnetism (AF) inside the superconducting 
phase under high pressure and magnetic field. 1 ' 2 ' 

CeCoIns is even a more fascinating case than CeRhlns 
as it was established by neutron 3 ) and NMR 4,5 ' exper- 
iments that for a field in the basal plane an AF phase 
referred as the Q-phase appears in the high field low 
temperature (HFLT) boundary pegged to H C 2(0), al- 
though above H C 2(0) the ground state is paramagnetic 
(PM). Prior to this recent experimental evidence, it was 
suggested 6-8 - 1 that the Fulde-Ferrel Larkin-Ovchinnikov 
(FFLO) state 9 ' 10 ' occurs in the HFLT region, with fi- 
nite momentum pairing created by the strength of the 
Pauli limit. Recent theoretical studies emphasize the key 
role of the Pauli depairing of the Cooper pairs 1 ^ and the 
possibility of a staggered field induced triplet component 
associated with the appearance of an AF phase coupled 
with or without a FFLO state 12 ' 13) 

Here we report on detailed magnetization M(T, H) ex- 
periments with the field H aligned along the c axis of the 
tetragonal crystal in order to search where are the singu- 
larities in the field dependence of the Sommerfeld coeffi- 
cient 7 of the specific heat and to determine if signatures 
of such a Q-phase are also present for this direction. For 
H along the c axis it was also suggested that H c2 (0) cor- 
responds to an AF quantum critical point. 14 ' 15 ) We have 
focused on the SC region close to H c2 (0) for H \\ c by 
using small field steps in order to look for evidence of 
the Q-phase. A HFLT phase has been reported to oc- 
cur along the H \\ c axis when H > 4.7 T and below 
T ~ 0.5 K based on NMR and fiSR experiments. 16 ' 17 ) 
However, no distinct transition line has been determined. 
The vortex lattice change from rhombic to hexagonal 18 ) 
close to H c2 {Q) at Hhex ~ 4.4 T is directly linked to the 
Pauli limit. 18 - 19 ) 
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Single crystals of CeCoIn5 were obtained by the In- 
flux method. Crystals grown previously in the same way 
in our laboratory have been used for NMR 5 ) and resis- 
tivity 20 ' measurements. 

The magnetization measurements were performed 
with a low temperature-high field SQUID magnetome- 
ter developed at Institut Neel. A miniature dilution re- 
frigerator allows for absolute value measurements by the 
extraction method at temperatures down to T = 70 mK. 
Notwithstanding, it is well known that large peak effects 
are present in this sample, and thus even a slight inho- 
mogeneity of the magnetic field, can give a large spurious 
signal when using the extraction method, as seen in a dif- 
ference between up and down extraction measurements. 
For our experimental setup, consistent results are ob- 
tained by averaging the up and down measurements to- 
gether. We digress to elucidate this subtle but important 
experimental detail: At H = 2 T, CeCoIns is near the 
maximum of the peak effect, and the field gradient in our 
magnetometer is approximately I Gauss over our normal 
scan length of 35 mm. At this field, we found that the 
difference between up and down extractions was approx- 
imately 0.007 /is/mole, or nearly the same magnitude 
as the averaged magnetization. (However, this difference 
becomes much smaller above 3.5 T where it is less than 
10% of the value of the magnetization, and drops to zero 
near 4.8 T.) During measurements of M vs T at fixed 
H, we have used two techniques in order to verify that 
averaging the up and down extractions is correct for this 
sample: One method is to reduce the scan length of the 
extraction. As the scan length is decreased, the difference 
between up and down extractions becomes smaller, and 
for scans below 15 mm, we could eliminate completely 
the difference, with the results converging directly on top 
of the averaged normal scan length data, although with 
more scatter in the data points. A second method is to 
place the sample in one of the detection coils of the mag- 
netometer, and then to vary the temperature without 
moving the sample at all, and thus measure the relative 
difference in magnetization at fixed field as a function 
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of temperature. (This technique was used for previous 
studies of the pseudo-metamagnetism of CeRu2Si2. 21 * 1 ) 
Again, these results (when correctly offset) fall on the 
averaged up and down extraction data. This gives us con- 
fidence in our method, and all the data shown in Figs. 1,2, 
and 3 are the averaged up and down measurements. 

Figure 1 shows M vs. H up to 8 T taken at T = 4.2 K 
and at T — 100 mK, well above and below the super- 
conducting transition at T c — 2.3 K in zero field. Good 
agreement is found with previous published magnetiza- 
tion data, 22 - 1 although the hysteresis and the observed 
peak effects in this sample are significantly smaller than 
previously reported. This is certainly due to the differ- 
ence between the field sweep technique of Ref.22 vs. the 
extraction method used here. But in addition, a part 
may be due to different sample qualities and their prepa- 
ration. The value of -ff C 2(0) = 4.96 T is in very good 
agreement with other bulk measurements like the specific 
heat, 7 ^ thermal expansion 23 ) and NMR. 16 ' The hystere- 
sis at H C 2(0) is quite small, 50 x 10~ 4 T and the jump of 
the magnetization is AM rj 3.5 x 10~ 2 /j,b, from an initial 
value of 0.1 [is just below -ff C 2(0). If the SC transition 
would be second order down to T = 0, a straightfor- 
ward extrapolation of the magnetization data suggests 
that H c2 (0) would be near 7 T. 

Our main interest in M(T) at constant field presented 
in this letter is that it is related to the field dependence of 
the specific heat coefficient C/T via the Maxwell relation 

d(C/T) _ d 2 M 
dH ~ dT 2 ' 
The high accuracy of the measurements associated with a 
careful determination of the reversible component makes 
possible an excellent estimate of d(C/T)/dH . In addi- 
tion, because it is a field derivative method of the spe- 
cific heat, it is quite sensitive not only to phase transi- 
tions but even to field induced crossovers. The Maxwell 
relation allows to obtain the bulk macroscopic magnetic 
field evolution of the specific heat, of course without the 



separation between the vortex and the supercurrent con- 
tributions. 

If the magnetization, assumed to be of purely elec- 
tronic origin, follows a characteristic Fermi-liquid tem- 
perature dependence M(T) f3(H)T 2 , the coefficient 
/3(H) is directly linked to the field derivative of the 
Sommerfeld coefficient dj/dH of the linear specific heat 
7(H)- 

Figure 2 shows the temperature dependence of M(T) 
at constant field above and below H &(()). The first ob- 
servation is that for fields H < H C 2(0) the magneti- 
zation decreases on lowering the temperature while for 
H > H C 2(0), M(T) increases on cooling. Clearly, the 
sign change in dM(H)/dT marks a singularity at H C 2- 
Another important observation is illustrated in Fig. 3 
where M is plotted vs. T 2 : for H < H C 2(0) and for low 
enough temperatures, M(T) is seen to follow a T 2 de- 
pendence. However, above H C 2(0) a dominant linear tem- 
perature dependence of M(T) has been observed down 
to the lowest measured temperature T = 70 mK. This 
last observation is still consistent with previous experi- 
ments where a Fermi-liquid regime was abserved above 
H C 2(0), but only at very low temperature T < 40 mK. 20 ' 
The clear observation of a T 2 dependence of M(H) for 
H < H C 2(0) demonstrates that the field dependence of 
the specific heat is governed by the growth of the linear 
term j(H) with magnetic field and not by other contri- 
butions. Similar effects have also been reported recently 
for H <C H C 2(0) when H is applied in the basal plane. 24 -* 
The huge field increase of 7 with field is associated with 
the gapless tendency of the SC on approaching H c2 (see 
Refs. 25-27) and the large values of 7 of the heavy fermion 
quasi-particles in their normal phase close to the mag- 
netic quantum critical point. 

From the T 2 law of M(T) measured below H c2 (0), 
/3(H) is determined, then dj/dH (see Fig. 4) and by 
integration the relative field variation Aj(H C 2 — e) ~ 
1.1 Jmole^K- 2 up to H c2 - e = 4.875 T. Aj(H) is 
shown in the inset Fig. 5(a). In the main panel Fig. 5(b), 
the data have been offset by a constant, adjusted in com- 
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Fig. 3. (Color online) Magnetization M as function of T 2 for dif- 
ferent magnetic fields H < H c 2(0) applied parallel to the c axis. 
Dashed lines indicate the a fit with a j3(H)T 2 dependence. 
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Fig. 4. (Color online) Field dependence ofdf/dH below H c 2(0). 
The arrow indicates the position of the vortex lattice phase tran- 
sition 18 ' at H hex = 4.4 T. In the insert we indicate H a « 4.1 T 
and Hp « 4.8 derived in Refs. 20,23. 



parison to direct specific heat data 28 ) at H — 3 T, giv- 
ing j(H C 2) ~ 1.2 Jmolc _1 K -2 . The field dependence of 
87 /8H does not show any distinct features that could 
be construed as the manifestation of a phase transition 
associated with the entrance into the Q state as observed 
for example in the field dependence of C/T for H || a at 
250 mK in Ref. 24. The absence of a Q-phase in CeCoIns 
for H || c is in good agreement with recent neutron scat- 
tering experiments claiming that a misalignment of 17° 
from the a axis to c axis leads to a collapse of the signal of 
the Q-phase. 29 - 1 However, we point out that a shoulder- 
like anomaly in the field dependence of 87 /8H occurs 
in the vicinity of the vortex lattice phase transition at 
Hhex ~ 4.4 T. In addition, different experimental tech- 
niques have indicated other characteristic critical fields 
associated with the interplay between SC and AF. Thus, 
for example, the location of the respective fields H a and 
Hp recently derived by thermal expansion and resistivity 




Fig. 5. (Color online) Inset (a): Relative field dependence of 
A-y(H) in the SC phase derived from magnetization measure- 
ments as explained in the text. Main frame (b): Comparison 
of 7(H) vs.H determined by different experiments. Solid line is 
7(H)vs.H at T = 250 mK from Ref. 30, solid diamonds from 
Ref. 15 for H > if c2 (0) taken at T = 120 mK and open di- 
amond for an extrapolation to 40 mK, for H < H c 2(0) taken 
from Ref. 28. Open circles indicate 7(H) from the present exper- 
iments and the adjustment of A7 has been made at H = 3 T. 
Triangles: 7 e (H) derived from resistivity experiments. 20 ) 



experiments are shown in the inset of Fig. 4. 20 ' 23 ) 

Above H C 2 the quasilinear T dependence of M(H) im- 
plies that the field dependence of C/T is very weak. This 
conclusion is in good agreement with the weak field de- 
pendence of C/T measured at T = 250 mK 30 ' but does 
not agree with the field dependence of C/T reported at 
T ~ 120 mK after the subtraction of the very large hy- 
perfine nuclear specific heat coming from Co and In nu- 
clei. 7, 31 > 32 ) Both elements Co and In have isotopes with 
large dipolar and quadrupolar moments. Thus a nuclear 
contribution has to be taken into account for the specific 
heat and magnetization analysis. Assuming a sole nuclear 
Zeeman term of the form AH 2 /T 2 , this additional term 
to 87 /dH gives a linear field correction 2AH/T 2 . The 
relative importance of the nuclear dependence will de- 
pend on the magnitude of the other contributions to the 
87 /dH derivative. Above H C 2 these latter contributions 
are small at least above 100 mK by comparison to its 
value at H C 2 — e, and the hyperfine contribution gives a 
significant channel for an increase in the T dependence of 
M on cooling leading to an apparent large temperature 
window where a linear temperature dependence of M(T) 
is observed. Thus we cannot derive just above H C 2(0) the 
variation of "f{H). To estimate this value in the normal 
phase at i? C 2(0), attempts have been made through the 
extrapolation of C/T to T = or through the indirect 
evaluation via the inelastic contribution to electron scat- 
tering. 

The precise determination of the A(H) coefficient of 
the T 2 term of the resistivity 20 - 1 (obeyed only below 
40 mK) for H > H C 2(0) confirms the quasi- divergence 
of A(H) close to i? C 2(0). 14 - ) In the physics of heavy 
fermion compounds the so-called Kadowaki- Woods re- 
lation A oc 7 2 is often obeyed due to strong local fluc- 
tuations. Scaling of \f A(H) with 7(H) predicts a value 
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of je(H C 2 +e) ~ 2.25 Jmole'^" 2 roughly in agreement 
with the low temperature extrapolation of the electronic 
contribution of C/T at 40 mK (see Fig. 5). A disconti- 
nuity of the 7 term at H c2 seems to occur in association 
with the first order nature of the metamagnetic phenom- 
ena at H c2 (0). Unfortunately, dHvA signals taken while 
varying the field through H c2 detect only the "light" 
branch. 33 ) However, the striking result is that the effec- 
tive mass increases from 15mo in the field range from 4.85 
to 4.95 T to 18m for a scan from 5 T to 5.2 T. A weak 
field dependence of j(H) appears for H > 7 T, the hy- 
pothetical value for a second order SC phase transition. 
Evidences that for H > 7 T another regime is achieved 
are also found in high field specific heat data. 34 ) 

The present results, for H < H c2 , clearly demon- 
strated that ■y(H) reaches its singularity at H c2 (0). There 
is no evidence of a Q-phase at least characterized by a 
singularity in the field dependence of 7 below if c2 (0), at 
least down to (H c2 - H)/H c2 > 0.02. 
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